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Abstract 

On the basis of the so-called "yukawaon" model, we found out a special form of the 
neutrino mass matrix M v which gives reasonable predictions. The M v is given by a multi- 
plication form made of charged lepton mass matrix M e and up-quark mass matrix M u . This 
M v has no adjustable parameters except for those in M e and M u . Here, M e and M u are de- 
scribed by one parameter a e (real) and two parameters a u (complex), respectively, and those 
parameters are constrained by their observed mass ratios. With this form of M v , in spite 
of having only three parameters, the M v can give reasonable predictions sin 2 29 atm c± 0.99, 
sin 2 26*13 — 0.015, Am| 1 /Am§ 2 — 0.030, (m ee ) ~ 0.0039 eV, and so on, by using observed 
values of mj/ifij,, m^/m T , m c /m t , and sin 2 9 so i ar as input values. 

PACS numbers: 14.60.Pq, 11.30.Hv, 12.60.-i, 

1 Introduction 

The observed masses and mixings of the quarks and leptons will provide a promising clue to 
a unified understanding of those fundamental particles. For such the purpose, not only inves- 
tigating a theoretical model, but also searching for phenomenological mass matrix seem to be 
still effective. As one of phenomenological mass matrix models of the quarks and leptons, the 
so-called "yukawaon" model [1] (a kind of "fiavon" models [2]) has been proposed. Here, the 
"effective" Yukawa coupling constants Yf are given by 

(Yf}) l3 = y J-((Y f ) l3 ) = 1,2,3), (1.1) 

A is an energy scale of the effective theory, and (Yf) are vacuum expectation value (VEV) 
matrices of scalar fields Yf with 3x3 components. (Hereafter we call fields Yf "yukawaons".) 
The most characteristic point in the yukawaon model is that all VEVs (Yf) are described in 
terms of only one fundamental VEV matrix (<& e ) oc diag( y/m e , •y/m^, ^/m T ). For example, in a 
yukawaon model with an 0(3) family symmetry [31 H], the charged lepton, neutrino, up-quark, 
and down-quark mass matrices M e , M u , M u and are given by 



M e Oc(y e )«($ e )($ e ), 

M v = m D M^m T D , 
m D oc (Y e ) , M R oc (Y r ) + {Y e ){Y e ), 

{Y R ) oc (<s> u )(p u ){Y e ) + (y e >(p u )($ u > + &,(($„) (p u >(y e ) + (y e )(p u )($ u )) 

M n oc (y u ) cx ($„)<$«), 

(* u ) OC ($ e )((P) +O u (X))(*e) ) 

M d cx (y d ) oc ($ e )((P> + a d (X))($ e ), 



where 
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in the diagonal basis of M e , while (P u ) is given by a form 



(Pu)u oc 




(1.10) 



in the diagonal basis of M u . 

In this paper, we will find out a special form of the neutrino mass matrix which is compatible 
with the observed neutrino data in spite of having no adjustable parameters. The form will be 
obtained along the lines of the yukawaon model by changing the structure of (Y e ) from Eq.(1.2). 

The neutrino mass matrix is still given by the form of seesaw type, M v = mDM^m^ with 
iri£) oc (y e ), but the Majorana neutrino mass matrix Mr is changed into a simple form 



Mr OC (Yr) OC {* u ){Y e ) + (Y e ){$u). 

Here, we assume that the charged lepton mass matrix M e is given by 

(Y e )<x(* )((E') + a e (X 2 ))(* ), 
differently from Eq.(1.2), and we also redefine as 

(* u )oc(* )((^) + OuW)(*o), 

where (E 1 ) = (E) = v E l and 



(1.11) 



(1.12) 



(1.13) 



(X 2 



vxS 2 = -vx 




(X 3 ) = v x S 3 = ^v x 




(1.14) 



Note that there are no term, no m 0jy term and no (P u ) in Eq.(l.ll), i.e. M v is simply 
given by 



In other words, there is no adjustable parameter in the present neutrino mass matrix, except for 
a e in M e and a u in M^ 2 (i.e. The purpose of the present paper is not to derive the mass 

matrix forms (1.11) - (1-13) theoretically, but to demonstrate that the phenomenological neutrino 
mass matrix (1.15) with Eqs.(l.ll) - (1.13) can be compatible with the present neutrino data in 
spite of quite few parameters. We predict sin 2 29 atm — 0.99, sin 2 2^13 ~ 0.015, Am^/Am^ — 
0.030, (m ee ) ~ 0.0039 eV, and so on, by using observed values of rrie/m^, m^/mr, and sin 2 9 so i ar 
as input values. In this paper, for simplicity, we do not discuss the down-quark mass matrix M d 
and the Cabibbo-Maskawa- Kobayashi [5] (CKM) mixing. 

In the next section, superpotentials for yukawaons and the assignments of the fields in the 
present U(3) yukawaon model are investigated. In Sec. 3, numerical results of the model are 
discussed. Sec. 4 is devoted to the concluding remarks. In Appendix A, ii-charge assignments 
are discussed. In Appendix B, we present a rotation matrix which transforms (A3} into (A2). 

2 Super potential 

In this section, we give superpotentials for the yukawaons and the assignments of the fields in 
the present U(3) yukawaon model. 

In the yukawaon model, the order of the fields is important. Therefore, in this paper, let us 
assume a U(3) family symmetry instead of 0(3) and denote fields 6* and 6 of U(3) as A and A, 
respectively. (Therefore, it should be noted that a term ABC is allowed, but ACB and BAC 
are forbidden.) In the U(3) model, for example, the relation (1.6) is re-expressed as 



with {E u ) = ugl. In order to distinguish each yukawaon from other yukawaons, although we 
assumed U(l)x charge in the 0(3) model [31 2], in this U(3) model, we assume only R charge 
conservation instead of U(l)x charge conservation. For the right handed neutrino sector (Yr), 
it should be noted that we cannot add Y e Y e term to Yr as in Eq.(1.4). In the old model, 
we assigned the U(l)x charges Qx only for gauge singlet fields, e.g. Qx{£) = Qx(H d ), i.e. 
Qx(Y e ) = —Q(e c ). Besides, we assumed Qx(e c ) = Q{v c ) in order to build a model without 
Y v . Therefore, we could obtain Qx(Yr) = Qx{Y e Y e ) in the old model. However, in this U(3) 
model, we cannot obtain R(Yr) = R(Y e EY e ). Besides, we cannot introduce a £„ term such as 
in Eq.(1.5). 

We assume the following superpotential W = Wy + W e + Wd + W' u + W u + Wr + We'- 




(1.15) 



(2.1) 



Wy = ^hYS C jH d + ^£ iY ^^H u + Xrv^v] 



+ ^u?YV qj H u + <fdCY* j gj H d + » H H u H d , 



(2.2) 
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Table 1: Assignments of R charges, where rx> = ry e — 2ro, r x = ryd — 2ro 
and R(<&u) = R(Yd) = ryd- The values in the third raw denote R charge values 
in a special case under the assumptions (A. 11) and (A. 14). For more details, see 
Eqs.(A.l) - (A.20) in Appendix A. 



W e = // e Tr[y e 9 e ] + ^Tr[$ (£' + a e X 2 )<5> ® e }, (2.3) 

W d = » d Tr[Y d e d ] + ^Tr[$ (£ + a d e ia «X 3 )$ & d ], (2.4) 

W' u = ^Tr[l> u e u/ ] + ^Tr[$ (£ + a u e ia "X 3 )<b Q u '], (2.5) 

W u = ^Tr[y u e u ] + ^Tt[$ u E u $ u Q u } + ^Tv[(E u E u fE u Q u ], (2.6) 

Wr = fj,RTr[Y R @ R ] + ^Tr[($ u £ u F e + F e ^$„) 9 R ], (2.7) 



where, in Eq.(2.2), g and t are SU(2)l doublet fields, and f c (f = u,d,e,v) are SU(2)i singlet 
fields. The other fields in Eqs.(2.3)-(2.7) have quantum numbers defined in Table 1. 

In Eq.(2.6), the third term has been added since it has the same R charge as that of Y U Q U . 
(See Eq.(A.20) in Appendix A.) The Xo u term plays a role in shifting eigenvalues of the up-quark 
mass matrix M u by a constant value. Details of R charge assignments are given in Appendix A. 

For the field E u , we assume an additional field E u , and consider a superpotential with a 

form 

W E = X E Tv[E u E u Q 8+1 } + \' E Tr[E u E u ] Tr[9 8+1 ], (2.8) 



where ©s+i is a field 8 + 1 of U(3) with (6s+i) = 0. The superpotential We leads to (E U )(E U ) oc 
1. We assume that the form 

(E u ) oc (E u ) = */ E diag(l,l,l), (2-9) 

is given by a specific form of the solutions {E U ){E U ) oc 1. 

In this paper, we do not discuss a superpotential which gives the observed charged lepton 
mass spectrum. We only use the observed charged lepton mass values as input values in (Y e ) e . 

Under the assumption that all fields take (0) = 0, SUSY vacuum conditions lead to VEV 
relation^ from Eqs.(2.3)-(2.7). That is, instead of Eqs.(1.2) - (1.8) in the previous model, we 
obtain the following mass matrix relations: 

M e oc (Y e ) oc ($ )(1 + a e S 2 )($ ), (2.10) 

M d oc (Y d ) oc ($ )(1 + a d e ia *S 3 )($ ), (2.11) 

(* u > oc <$ )(1 + a u e ia "S 3 )(^ ), (2.12) 

M u oc (Y u ) oc ($ u ) • 1 • ($ tt ) + (^ U ) 2 C„1, (2.13) 

M fl oc (Yr) oc (* u ) • 1 • (Y e ) + (Y e ) • 1 • ($ tt ). (2.14) 
Here, we can take a diagonal basis of ($o) without loosing the generality: 

($o) = diag(«i,«2,«3) = -vodiag(xi,x 2 ,2;3), (2.15) 

where we have normalized + x| + x 2 = 1. The neutrino mass matrix M v is given by 

a seesaw type M v = moMZ rn^ with tub oc M e similar to Eq.(1.4) [but there is no m^]. 
Note that the previous relations (1.2) - (1.8) were given at a diagonal basis of the VEV ($ e )> 
while present relations (2.10) - (2.14) are given at a diagonal basis of the VEV ($o)- Here the 
numerical matrices S3 and S2 are defined by 



53 = 1 



/I 1 1 
1 1 1 

\1 1 1 



1 1 1 I , S 2 = - 



/ 1 1 

1 1 I , (2.16) 
V 



at the diagonal basis of the VEV (^o)- Note that the VEV matrix (Y e ) in Eq.(2.10) is no more 
diagonal in this basis. 

In obtaining the mixing matrices, the common coefficients are not important. Here we 
have taken ve< = VX2 and ve = vx3 for simplicity. The £ u term in Eq.(2.13) comes from the 
new term given in Eq.(A.20). This term contributes to the up-quark mass ratios, while not to 
the up-quark mixing matrix, so that it does not change the predictions for the neutrino mixing 
parameters. We suppose that the contribution from such the higher dimensional term (A. 20) 
is considerably small, so that it also does not visibly affect the up-quark mass ratio m c /mt, 
although it can slightly affect m u /m c . 



1 For example, in obtaining the relation (2.10), we have assumed a vacuum with (Q e ) = 0, so that the 
conditions 8W/dY e = and dW/d&o = do not affect other VEV relations obtained from SUSY vacuum 
conditions OW/OOa = (A/e). We assume that the observed SUSY symmetry breaking is induced by a gauge 
mediation mechanism (not including family symmetry), so that our VEV relations among yukawaons are still 
valid in the quark and lepton sectors after the SUSY is broken . 



3 Numerical results in the up-quark and neutrino mass matrices 



In this section, we investigate whether the new VEV matrix relations (2.10) - (2.14) can well 
describe the observed neutrino mixing parameters together with the observed up-quark mass 
ratios or not. 

Since the charged lepton mass matrix given by Eq.(2.10) is not diagonal, the lepton mix- 
ing matrix [Pontecorvo-Maki-Nakagawa-Sakata (PMNS) [7] mixing matrix] U in the present 
conventions is defined by 

U = ul L U uL , (3.1) 

where U e i and U u l are defined by 

U? L (Y e )U eL = (Y e ) dmg = (Y e ) e , (3.2) 



ul L {MtM v )U vL = (MtM u ) dia \ (3.3) 



and M v is given by 



M - = f|(^r) omr 1 ^). (3.4) 

Neutrino mixing parameters we discuss are tan 2 9 so i ar = \U\2\ 2 /\Un\ 2 , sin 2 26 atm = 4| 1 2 1 C^33 1 2 , 
and |?7i3| 2 . Here Uij are the matrix elements of the lepton mixing matrix defined by (3.1). 
The matrix M u in (2.13) is diagonalized as 

ul L (MlM u )U uL = (MlM u f a °. (3.5) 

Here U u l is a mixing matrix among left-handed up-quarks uu ■ (In the present paper, the mass 
matrices (i.e. (Yf)) are defined by Eq.(2.2). Therefore, the conventions of the mixing matrices 
are somewhat changed from the conventional ones.) Note that since the VEV matrix is 
complex and (Y u ) is given by Eq.(2.13), the diagonalization of the up-quark mass matrix must 
be done by Eq.(3.5). 

3.1 Parameters in the model 

The mass matrices for quarks and neutrinos in the 0(3) model have been described in 
terms of the fundamental VEV matrix (3> e }. On the other hand, the fundamental VEV matrix 
in the present model is ($0) defined by Eq.(2.10) in which we have new parameter a e . Thus the 
number of parameters are increased by one compared with the previous model (1.2). On the 
other hand, we cannot bring neither the £ u term given in Eq.(1.5) nor (P u ) u defined in Eq.(l.lO) 
into the present model, so that there are no parameters which are corresponding to and P u . 

The VEV of (<J»o) = diag(i>i, 112, V3) is related to the charged lepton mass matrix M e as 
follows: 



( (l + \a e ) v\ \a e viv 2 

\a e viv 2 (l + \a e ) v\ (3.6) 



where k = — (\ e / jj, e A)(y e (H a l) / A)vx' , so that we obtain 



m e + = k \1 + -a e j (vf + v%), (3.7) 

iTieiri/j, = k 2 (l + a e )vfv 2 , (3-8) 

and m T = fet>|. Here, since we are interested only in the relative ratios among the eigenvalues of 
the charged lepton mass matrix M e , the common coefficient k is not a parameter of the model. 
The 3 parameters a e , v\/v2 and v-ijvj, are sufficient to determine the two charged lepton mass 
ratios and charged lepton mixing matrix U e which is described only by one parameter 9f 2 - [The 
mixing angle 9\ 2 is not observable. The observed quantities are parameters of the lepton mixing 
matrix defined by Eq.(3.1).] Therefore, when we give a value of the parameter a e , the values of 
Vi are completely determined by the input values of the charged lepton masses. In other words, 
even when we give three charged lepton masses as the inputs, one of the free parameters still 
remains. 

Thus, in the present model, we have 4 parameters a e , a u , a u and Cm (except for the input 
values m e , m^, and m T ) for the up-quark and neutrino mass matrices. On the other hand, the 
number of the predictable quantities are 12, i.e., 2+2+2 mass ratios (up-quark, charged lepton 
and neutrino mass ratios) and 4+2 PMNS mixing parameters (including two Majorana phases). 
At present, we know 6 observed values of y/m u /m c , y/m c /mt, tan 2 6 so i ar , sin 2 2# atm , sin 2 2#i3, 
and R u = Am^/ Am 2 2 in addition to the charged lepton masses. 

The term Q u \ in M u given in Eq.(2.13) does not affect the up-quark mixing matrix, so 
that it also affects neither quark or lepton mixing matrices. Since we suppose \Cu\ 2 <C 1, the 
term almost does not affect y/m c /mt, although it can slightly affect \fm u /m c . As a result, 
the present model predicts 11 observables by using the three parameters (a e ,a u ,a u ). In other 
words, the value of \J m u /m c is not " prediction", and it is a quantity which can be adjustable 
by the additional parameter Q u freely. 

3.2 Numerical results 

Now let us show the results of numerical analysis of the model. First, we show, in Fig. 1, 
the a u dependences of the quantities y/m c /m t , tan 2 9 so i ar , and sin 2 26 atm with taking typical 
values of a e = 3, 30, 100 and a u = 0°, 15° in order to see rough parameter behaviors. 

As seen in Fig. 1, we can find that (i) the value of \fmjm~t takes a maximum value at 
a u ~ —3 insensitively to the values of a e and a u ; (ii) since the maximum value of sin 2 29 a t m 
shows sin 2 26 a t m ~ 1 which is in favor of the observed value, we must search for a parameter 
set (a e , a u ,a u ) which gives a maximum value of sin 2 29 a t m ; (iii) a case with a small value of a e 
gives a large value of tan 2 9 so i ar compared with the observed value tan 2 9 so [ ar ~ 0.5 (see Fig. 1 
(a)), so that such a case is ruled out; on the other hand, a case with a large value of a e gives a 
small tan 2 9 so i ar (see Fig. 1 (c)), so that such a case is also ruled out; (iv) as a result, a region 
of (a e , a u ) which can give sin 2 29 a tm — 1 and tan 2 9 so i ar ~ 0.5 is (a e , a u ) ~ (30, —3). 

Next, in order to determine parameter values (a e ,a u ,a u ), let us illustrate, in Fig. 2, the 
a u behaviors of m c /m t , tan 2 9 so i ar and sin 2 29 atrn at a e ~ 28 and a u ~ —3. From Fig. 2, we 




Figure 1: y/m c /mt, tan 2 9 so i ar , and sin 2 29 a t m versus a parameter a u for typical parameter values 
a e = 3 (Fig. 1 (a)), a e = 30 (Fig. 1 (b)), and a e = 100 (Fig. 1 (c)) with a u = 0° (solid curves) 
and a u = 15° (dashed curves). Curves "r23", "solar", and "atm" denote "r23"= \fm c /m t x 10, 
"solar" = tan 2 6> so ; ar , and "atm"= sin 2 2# a4m , respectively. 




(a) (b) (c) 



Figure 2: y/ m c /mt, tan 2 9 so i ar , and sin 2 29 a t m versus a parameter a u for typical parameter 
values a e = 26 (Fig. 2 (a)), a e = 28 (Fig. 2 (b)), and a e = 30 (Fig. 2 (c)) with a u = -2.9 
(dashed curves), a u = —3.0 (solid curves), and a u = —3.1 (dot-dashed curves). Curves "r23", 
"solar", and "atm" denote "r23"= y/m c /m t x 10, "solar" = tan 2 9 so i ar , and "atm" = sin 2 29 a tm , 
respectively. 

search for the value a u which gives the observed value [8] yjm c /mt = 0.0600^q'qq47 at fi = mz- 
We find that the value a u ~ 15° can give a reasonable fit y/m c /mt = 0.0600 insensitively to the 
other parameters. 

Therefore, by fixing the value a u = 15°, we illustrate the contour lines of ^m c /m t and 
tan 2 9 so i ar in the (a e ,a u ) plane in Fig. 3. The curves denote (a e ,a u ) which gives the observed 

values ^m c /m t = 0.0600±g;oo47 and tan2 e foiar = °- 47 -ao3 M- As seen in Fi §- 3 - we have two 
intersection points of the curves of \Jm c jmt and tan 2 so [ ar . For the center values y/m c /mt = 
0.060 and tan 2 9 so i ar = 0.47, the solutions (a e ,a u ,a u ) are 



(26.7, -2.88, 15°), (27.9, -3.18, 15°). 



(3.9) 




Figure 3: Contour lines of y/ m c /mt and tan 2 9 so iar in the (a e , a u ) plane in the case of a u = 15°. 
As input values, x/m c /m t = 0.06001^7 and tan 2 9° b J ar = 0.47±g;g| have been used. Curves 
with dot-dash, solid, and dash denote the upper, center, and lower observed values, respectively. 
Curves "r23" and "solar" denote "r23"= \/m c /m t x 10 and "solar" = tan 2 9 so i ar , respectively. 



Input (a e ,a u ,a u ) 


\Jm c /mt 


tan 2 9 so l ar 


sin 2 29 atm 


sin 2 26»i3 


R v 


(26.7,-2.88,15°) 


0.0603 


0.470 


0.990 


0.015 


0.0303 


(27.9,-3.18,15°) 


0.0601 


0.469 


0.980 


0.011 


0.0204 


upper 


+0.0045 


+0.05 






+0.0027 


Observed value 


0.0600 


0.47 


> 0.92 [9J 


< 0.15 [9] 


0.0312 


lower 


-0.0047 


-0.03 






-0.0023 



Table 2: Predicted values for the parameter values (a e ,a u ,a u ). 



We list our prediction values for these parameter solutions in Table 1. Of the two solutions 
obtained from the input data \J m c /mt and tan 2 9 so i ar , Table 2 suggests that we should take the 
former one considering the observed value of R v [5] 



Am2 solar 
Am ltm 



(3.12 



+0.27N, 
-0.23 > 



X 10 



-2 



(3.10) 



For reference, we also illustrate the behavior of predicted values for input values (tt e ,(x w ,cy w ) 
around the parameter solutions (3.9) in Fig. 4. As seen in Fig. 4, the predicted values sin 2 29 a t m 
and tan 2 9 so i ar are insensitive to the parameter values a e , a u , and a u around the values (a e , a u ,a u ) = 
(26.7,-2.88,15°). However, y/ m c /mt and |C/i3| 2 (and also i?„) are somewhat dependent on 
these parameters. Since these parameter values are mainly obtained by taking the input value 
\J m c /m t = 0.0600, if the input value changes, then the predicted values will also change. 

So far, we have not discussed the value of m u /m c . In the present model, the value of 
"m u /m c is always adjustable by the parameter £ u given in Eq.(2.13) without affecting other 




Figure 4: Predicted values versus a parameter a u for a e = 26.7 and a u = —2.88. Curves 
"r23", "solar", "atm", "ul3", and "R" denote "r23"= yjm c /m t x 10, "solar"= tan 2 e solar , 
"atm"= sin 2 26> afm , "ul3"= |£7i 3 | 2 x 100, and "R"= R v x 10, respectively. For reference, curves 
for {a e ,a u ) = (26.7, -3.00) (dash curve), (26.7, -2.80) (dot curve), (29, -2.88) (dot dash curve), 



and (25,-2.8 


8) (2-dot dash curve) are 


illustrated 


in addition to the curve (solid) for (26.7, — 2.81 




Input (u 




tan 2 9 so l ar 


sin 2 29 atm sin 2 26>i 3 5 C p ( J) 


R u 




3.8 x 10~ 7 


0.0180 0.0603 
0.0454 0.0602 


0.470 
0.470 


0.990 0.015 -102.5° (-0.0139) 
0.990 0.015 -102.5° (-0.0139) 


0.0303 
0.0303 



Table 3: Predicted values versus ( u parameter. Other parameters (a e ,a M ,a M ) have 
taken the same values (26.7, —2.88, 15°) as those in Table 2. 



predicted values. In order to fit the predicted value of \J m u /m c to the observed value [8] 
y/m u /m c = 0.0453^q'q];q, we choose ( u as ( u = 3.8 x 10~ 7 . As seen in Table 3, the value of ( u 
almost does not change the numerical predictions given in Table 2. 
In conclusion, we take the parameter set 

(ae, a u , Ou, Cu) = (26-7, -2.88, 15°, 3.8 x 10^ 7 ). (3.11) 

Then, we predict neutrino masses 

m„i ~ 0.012 eV, m u2 ~ 0.015 eV, m u3 ~ 0.051 eV, (3.12) 

by using the input value Am| 2 — 0.0024 eV 2 . We also predict the effective Majorana mass (m ee ) 
in the neutrinoless double beta decay [TO] 

(mee) = \miU^ + m 2 C/ e 2 2 + m 3 U% 3 \ ~ 0.0039 eV. (3.13) 

It is worthwhile noticing approximately degenerate neutrino masses m v \ ~ m U 2- 



4 Concluding remarks 



In this paper we have found out a special form of the neutrino mass matrix M v based on a 
yukawaon model with U(3) family symmetry, which has quite few free parameters. With this 
form of M v , the M v can give reasonable predictions in spite of having no adjustable parameters, 
i.e. the M v is simply given by the form (1.15), and the mass matrices M e and M u include 
parameters a e and a u , respectively, which are fixed by their observed mass ratios. In this 
yukawaon model, the yukawaon VEV matrices are described in terms of a new fundamental 
VEV matrix ($o)- For example, the yukawaon VEV matrix (Y e ) for the charged leptons is given 
by (2.10) which has the structure of (1 + a e S2) with a new parameter a e . This structure in (Y e ) 
has been chosen from a phenomenological point of view and there is no reason why (Y e ) takes 
such a form. Nevertheless if we accept the form (2.10), then we can obtain a simple form of VEV 
matrix (Yr) f° r right-handed neutrinos without introducing the somewhat strange VEV matrix 
P u and term that were introduced in the 0(3) model to get the observed nearly tribimaximal 
neutrino mixing [llj . 

The new model has only four parameters (a e , a u , a u , £ u ) as far as the up-quark and lepton 
sectors are concerned, on the other hand, we have 12 observable quantities (2 up-quark mass 
ratios, 4 lepton mass ratios, and 4+2 lepton mixing parameters). The parameter Q u affects only 
the prediction of m u /m c , so that we have fixed it by the observed value of r" 2 = y/m u /m c . 
The parameter a u is sensitive only to m c /mt, so that we have fixed by the observed value of 
r 23 = \l m c/ m t as seen in Fig. 2 (b). The parameter a e is determined from the cross point of 
the predicted values of and tan 2 9 so i ar in the a e -a u plane. Note that we have used only the 
observed values of r% 3 and tan 2 9 so i ar in order to fix the three parameters (a e , a u , a u ). Although 
we have tacitly used sin 2 26 atm ~ 1 , we have not used the observed value of sin 2 29 a tm explicitly. 

On the other hand, for the remaining 2 down-quark mass ratios and 4 CKM mixing param- 
eters, we have additional 2 parameters (o^ and a^). Regrettably, we cannot obtain reasonable 
predictions with the two parameters, although we can fit the values of down-quark mass ratios 
and V us . The situation is the same as in the previous 0(3) model. We must introduce a phase 
matrix with two parameters (4>i,4>2) and a common mass shift term modi. Then, five pa- 
rameters can fit six observables barely. Therefore, the model is not so attractive for down-quark 
sector. In this paper, we did not demonstrate the explicit numerical fitting for down-quark mass 
rations and CKM mixing parameters. 

The present U(3) model have the following interesting features in the lepton sector: 
(i) The model predicts sin 2 28±^ ~ 0.015. In the previous 0(3) model, the predicted value 
of |^13| 2 w& s invisibly small, i.e. | tT" 13 1 2 ~ 10~ 4 . The T2K experiment [12] put a constraint 
0.03 < sin 2 26<i3 < 0.28 (90% C.L.) for 5cp = and a normal hierarchy. Our predicted value 
sin 2 2#i3 = 0.015 seems to be somewhat lower than the experimental lower bound. However, as 
seen in Table 3, our prediction on 5c p gives 5c p = —103°, which decreases the lower bound 
0.03 of the T2K result to 0.02. Besides, the Double CHOOZ experiment |13j has reported that 
sin 2 20i3 = 0.085 ± 0.029 ± 0.042 at 68% CL. The lower value is sin 2 2% = 0.014. Therefore, 
we consider that the predicted value sin 2 2#i3 = 0.015 is yet not ruled out, although the status 
is considerably severe. 



(ii) It also predicts a reasonable value of R v = Am^/Am^ ~ 0.03 in contrast to the case 
of the 0(3) model in which we could not predict R u . (In the previous model, the value of R v 
needed to adjust the additional free parameter m^ 1 in Eq.(1.4).) 

(iii) The present model gives approximately degenerate neutrino masses m v \ ~ m v %. The 
predicted value for the effective Majorana mass (m ee ) — 0.0039 eV in the neutrinoless double 
beta decay will be within our reach of the future experiments. 

The big ansatz is the existence of the X2 term in the charged lepton sector (1.12). At 
present, there is no idea on this term. Besides, it seems that the present lepton mass structure 
is ill matched with the charged lepton mass relation [H] 



The purpose in the early stage of the yukawaon model was to predict the charged lepton mass 
relation (4.1). The bilinear form (1.2) for the charged lepton mass matrix was indispensable to 
predict |15| the relation (4.1). If we adopt the present scenario, we must reconsider the origin 
of the charged lepton mass spectrum. However, in this paper, we do not use the relation (4.1), 
but only use the observed charged lepton mass values as input values. Therefore the bilinear 
form such as (1.2) is not necessarily required in this paper. Nevertheless, the formula (4.1) 
is still attractive. On the other hand, it is also attractive that we can predict 12 observables 
(2+2+2 lepton and up-quark mass ratios, and 4+2 PMNS mixing parameters) under 4 adjustable 
parameters (a e , a u , a u , ( u ) if once we accept this ansatz (1.12). It is a future task how to 
understand the existence of X2 term. 

In conclusion, although the form M u is, at present, not one which is derived from a rigid 
theoretical ground, the form will offer a suggestive hint for a unification model of quark and 
lepton mass matrices. 
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Appendix A: R charge assignments 



In the present model, as well as in the 0(3) model, we construct a model without introducing 
a yukawaon Y v by replacing Y v by Y e . The simple way to guarantee that the yukawaon Y e couples 
not only to the charged lepton sector but also to the Dirac neutrino sector is to introduce the 
following R charge assignment, 

R{u c ) = R(e c ) = r e , (A.l) 

R{H U ) = R{H d ) = 1. (A.2) 
The R charge of (E U E U ) is free parameter in the form (2.8). For simplicity, we take 

R(E U E U ) = 22(08+1) = 1. (A3) 

Hereafter, we will denote R(E U ) and R{E U ) as te and 1 — ?e, respectively. Each yukawaon is 
distinguished from other yukawaons by the R charges. If we define a parameter n as 

n = 2[R{Y R )-R(Y e )l (A4) 

then, we can express the R charges of the other fields from Eq.(2.2) as follows: 

R{() =r e + \{n - 2), (A5) 

i?(y e ) = I(4-n)-2r e , (A6) 

R(Y R ) = 2 - 2r e , (A7) 

R(Y u ) = n-l + r E , (A8) 

R(Y d ) = ^(n-2) + fE, (A9) 

R(u c ) + R(q) = 2 - n - r E , (AlO) 

R(d c ) + R(q) = 2 - l -n - f E . (All) 
From Eqs.(2.6) and (2.7), we obtain 

R(Y U ) = 2R($ U ) + R(E U ), (A12) 

R(Y R ) = R($ u ) + R(Y e ) + R(E U ), (A13) 
respectively. From Eqs.(A.12) and (A. 13), we obtain a relation 

R(Y U ) = n- R{E U ) =n-l + f E . (AAA) 

The relation (A. 14) leads to 

R(Y U Q U ) = (n - l)R{E u E u ) + R(E U Q U ). (A15) 



Only when the value n is a positive integer, Eq.(A.15) means that an additional term 



j^^Tr[(E u E u ) n ~ 1 E@ u ], (A.16) 

can appear in the expression (2.6). Note that if n is not a positive integer, the factor (E u E u ) n ~ 1 
does not have a physical meaning, because a term with (E u )~' 1 cannot appear in the superpo- 
tential terms. Therefore, the n defined in Eq.(A.4) is allowed only for n = 1, 2, • • • . 

As we see in Eqs.(A.5) - (A. 11), these R charges are described by four parameters r e , R(q), 
Te and n. Therefore, in order to fix these R charge values, we have to assume four constraints 
for these R charges. On the other hand, the fields Y e , Yr, Y u , Y^, $ n , and E u are gauge singlets, 
so that they must be distinguished only by R charges. We can choose a suitable parameter set 
(n,r e ,r q ,fE)- Here, let us demonstrate an example of R charge assignments, although it is not 
the purpose of the present paper to give such an explicit R charge assignment. 

For example, we put the following working hypothesis: 

R(Y V ) + R(Y e ) = 0, R(Y U ) + R(Y d ) = 0, (A 17) 

R(u c ) + R(d c ) = 0, R(u c ) + R(e c ) = 0. (A18) 

The constraint (A. 17) is an analogy that the Yukawa coupling constants in the standard model do 
not have R charges. The constraints (A. 17) and (A. 18) leads to the relation R(£) = R(q) = 1- Of 
course, since the yukawaon Y v has been replaced by Y e in the present model, the first constraint 
in Eq.(A.17) reads as R(Y e ) = 0, and since R{v c ) = R(e c ) in the model, the second constraint 
in Eq.(A.6) reads as R(e c ) = 0. Since R(Y e ) is given by Eq.(A.7), the requirement R(Y e ) = 
together with R(e c ) = requires n = 4. Thus, the constraints (A. 17) and (A. 18) fix the 
parameters (n,r e ,r q ,fE) as 

n = 4, R(e c ) = 0, R(q) = 1, R(E U ) = -2. (A19) 

The explicit values of these R values are listed in Table 1. Since the R charges of X' and X 
are still free parameters, we take i?($o) = \ for simplicity. As we see in Table 1, the fields Y e , 
Yr, Y u , Yd and E u can safely have different R charges from each other. 

Thus, the assumption can lead to plausible R charge values (A. 19), so that we consider 
that the assumption is reasonable. Now we have an additional term, 

^Tv[(E u E u fE u @ u ], (A20) 
which should be included in M u given in Eq.(2.6). 



Appendix B: Rotation from S3 into S2 



We define 



.5', 




So 




(B.l) 



which are invariant under permutation symmetries S3 and S2, respectively. A rotation matrix 
R which transforms the matrix S3 into S2 has been discussed in Ref.|16j. The rotation matrix 
R is given by 

RS3R? = 5*2, 



where R is defined as follows: 



7T . 



R = R3(-j)TR 3 (0)A, 
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The matrix A is known as a matrix which diagonalizes the matrix S3 into 



AS,A T 
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and also 
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The explicit form of R is given by 



where s = sin 9 and c = cos 6. 

In the expression (B.9) of R, when we define 
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the rotation matrix R is expressed as follows: 

( 

R= yjl-zi 

\ ^/f(«3 - ^2 

Here, (i = 1,2, 3) satisfies 
and we can choose Z{ such as 

z\ + z 2 + 2:3 

Suggested from the charged lepton mass relation (4.1), if we choose z, 



^m el + m e2 + m e3 ' 



where (m e i,m e2 ,m e3 ) = (m e , m^, m T ), then, the matrix i? satisfies 
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where 



zs*z = I 
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the following relation holds: 

RZS 3 ZR T = -Z u (B.18) 
3 

where Z\ is defined by Eq.(B.8). However, note that, from Eqs.(B.8) and (B.18), we cannot 
conclude RZ = (l/y/3)TA. 

Thus, it seems the rotation matrix R from S3 into £2 is deeply related to the charged lepton 
mass relation (4.1), but it is not clear why the form S2 appears in the charge lepton sector. This 
is still an open question at present. 
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